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ABSTRACT 

Dust polarization orientations in molecular clouds often tend to be close to tan- 
gential to the Stokes / dust continuum emission contours. The magnetic field and the 
emission gradient orientations, therefore, show some correlation. A method is pro- 
posed, which - in the framework of ideal magneto-hydrodynamics (MHD) - connects 
the measured angle between magnetic field and emission gradient orientations to the 
total field strength. The approach is based on the assumption that a change in emis- 
sion intensity (gradient) is a measure for the resulting direction of motion in the MHD 
force equation. In particular, this new method leads to maps of position-dependent 
magnetic field strength estimates. When evaluating the field curvature and the gravity 
direction locally on a map, the method can be generalized to arbitrary cloud shapes. 
The technique is applied to high-resolution (~ O'.'l) Submillimeter Array polarization 
data of the collapsing core W51 e2. A tentative ~ 7.7 mG field strength is found when 
averaging over the entire core. The analysis further reveals some structures and an 
azimuthally averaged radial profile ~ r" 1 / 2 for the field strength. Maximum values 
close to the center are around 19 mG. The currently available observations lack higher 
resolution data to probe the innermost part of the core where the largest field strength 
is expected from the method. Application regime and limitations of the method are 
discussed. As a further important outcome of this technique, the local significance of 
the magnetic field force compared to the other forces can be quantified in a model- 
independent way, from measured angles only. Finally, the method can potentially 
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also be expanded and applied to other objects (besides molecular clouds) with mea- 
surements that reveal the field morphology, as e.g. Faraday rotation measurements in 
galaxies. 

Subject headings: ISM: clouds — ISM: magnetic fields, polarization — ISM: individ- 
ual (W51 e2) — Methods: polarization 



Introduction 



Mag netic fields are being rec ognized as one of the key components in star formation theo- 
ries (e.g., iMcKee & Ostrikerll2007|) . Nevertheless, their exact role in the formation and evolution 
of molecular clouds is still a matter of debate in the literature. On the observational side, one 
of the difficulties is to accurately measure the magnetic field strength. The Zeeman effect still 
provides the only known method o f directly measuring magnetic field strengths along the line 
of sight in a molecular cloud (e.g., iTroland & Crutcherj 120081: ICrutcher et al.ll2009|) . A statistical 
method, aimin g at separating a three-d imensional uniform field and a non-uniform field compo- 
nent is given in JokiDii & Parker (11969|) and was later extended and applied to optical polarization 



nent is given rnlJOKipn & rarkeri ana was later extended and applied to optical poianzatioi 

data in lMyers & Goodman! dl 99 ll). The Goldreich-Kylafis e ffect, the linear polarization in molec 



ular line spectral emission (|Goldreich & Kylafislll98lLll982|) . is a mechanism to trace the magnetic 
field morpholog y. Up to now, the field structures of on ly a few sources have been detected with 
this effect (e.g.. ICortes et al]|2008l : iBeuther et al.l 120 10f) . Polarization of dust thermal emission at 
infrared and submillimeter wave lengths provides another method to study magnetic field proper- 
ties (e.g. jHildebrand et al. 2000 ). The dust grains are thought to be aligned with their shorter axis 



parallel to the magne tic field lines. Therefore, the emitted light appears to be polarized perpendic - 
ular to the field lines dCudlip et al.ll982l : Iffildebrand etaD 19841 lHildebrandjl988[ iLazarianlboOol) . 



Radiative torques are likely to be responsible for the dust alignment ( Draine & Weingartnerl 1996i 
19971 : lLazarianll2000l : Icho & Lazarianll2005l : Lazarian & Hoangll2007l) . 



Dust polarization obs ervations have been obtained both f or cores together with their surround- 
ing diffuse r material (e.g., Schleuningl 1998[ Lai et al. 2001 ) and for higher-resolution collapsing 



cores (e.g., Girart et ahbood : Tang et al. 



2009J). Complementary to Zeeman splitting, dust polariza- 



tion measurements are probing the plane of sky projected magnetic field structure. They, however, 
do not directly yield a magnetic field strength, but only some field morphology. Additional model- 
ing is needed in order to derive the actual field strength perpendicular to the line of sight. A three- 
dimensional model of an hourglass magnetic field configuration is fitted to dust observations in, e.g. 
Kirbyl (|2009|) to yield a mean field strength. With pinched field configurations det ected, comparing 
gravitating force and field tension gives an upper limit for the field strength (e.g.. lSchleuningll998 : 
Tang et al.ll2009c iRao et al.ll2009|) . Another commonly used technique is the Chandrasekhar-Fermi 
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(CF) method (|Chandrasekhar & Fermill 1953b . where a field strength is derived from the field dis- 



persion (around a mean field direction) combined with a turbulent velocity infor mation. Various 



Houde 2004: 


Falceta-Goncalves et al. 


2008 


) and tested through magneto-hydrodynamics (MHD) 


simulations ( 
The work by 


Ostriker et al. 


2001 


: Padoan et al. 2001: Heitsch et al. 


2001: Kudoh&Basu 


2003). 


Hildebrand et al. 


(2 


009) further refines this method - noting that a globally derived 



dispersion might not reflect the true contribution from magneto-hydrodynamic waves and/or turbu- 
lence - by introducing the concept of a dispersion function (structure function) around local mean 
magnetic field orientations. 

However, the CF method and all its variations are statistical in nature, relying on (large) 
ensembles of polarization segments. Thus, they typically provide one ensemble-averaged value 
for the magnetic field strength for an entire map. Consequently, zones of weaker and stronger 
field strengths can not be separated. It is, therefore, non-trivial to evaluate the relative, generally 
position-dependent role of the magnetic field (compared to, e.g., gravity, turbulence). The goal 
of this paper is to introduce a new method which leads to a local magnetic field strength and, 
therefore, can map the position-dependent field strength over an entire molecular cloud 

The paper is organized as follows. Section [2] summarizes the Submillimeter Array (SMA) 
polarization observation of the source W51 e2 which serves as an illustration case. Section [3] 
introduces our new method, starting with an observed correlation, the dust - MHD connection and 
a qualitative analysis of the predicted features of the method. A possible extension of the method 
to an arbitrary cloud shape is discussed in Section HI The results are presented in Section [5J We 
follow up with a discussion in Section [6] which addresses application regime and limitations of the 
method, and which also proposes a measure to evaluate the dynamical role of the magnetic field in 
different zones of a molecular cloud. Summary and conclusion are given in SectionfTJ 



2. An Illustration Case: W51 e2 



The W5 1 e 2 map (Figur e [B analyzed in the following sections was obtained with the SMA 
interferometer jTang et al. 20091) . The observation was carried out at a wavelength of 0.87 mm 
where the polarization of dust thermal emission is traced. In the extended array configuration, 
the synthesized beam resolution was about O'.'l with a primary beam (field of view) of about 30" 
at 345 GHz. A detailed descrip tion of the data analysis with a discussion and interpretation of 
the obtained images is given in iTang et al.l (|2009|) . The re-constructed features depend on the 
range of Mv-sampling and weighting. Nevertheless, these data currently provide the highest angular 
resolution (6>) information on the morphology of the magnetic field in the plane of sky obtained 
with the emitted polarized light in the W51 star formation site. This study is part of the program on 
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the SMALU (H o et al] |2004) to study the structure of the magnetic field at high spatial resolutions. 



W51 e2 is one of the strongest m m/submm continuum sources in the W51 region. Located a t 
a distance of 7 kpc (IGenzel et al.lll98l[) . 1" is equivalent to ~ 0.03 pc. IChrysostomou et al.l (I2002|) 
measured the polarization at 850/im with SCUBA on JCMT across the molecular cloud at a scale 
of 5 pc with a binned resolution 6 rs 9."3. At this scale, the pol arization s hows a morphology that 
changes from the dense cores to the surrounding diffuser gas. lLai et al.1 (1200 1|) reported a higher 
angular resolution (3") polarization map at 1.3 mm obtained with BIMA, which resolved out large 
scale structures. In contrast to the larger scale polarization map, the polarization appears to be 
uniform across the envelope at a scale of 0.5 pc, enclosing the sources e2 and e8. In the highest 
angular resolution map obtained with the SMA at 8 70 jum with 9 fa O '.'l ', the polarization patterns 
appear to be pinched in e2 and also possibly in e8 (|Tang et al.1 12009). In particular, the structure 
detected in the collapsing core e2 reveals hourglass-like features. A statistical analysis based on 
a polarization structure function (of second order) shows a turbulent to mean magnetic field ra- 
tio which decreas es from the larger (BIMA) to the smaller (SMA) scales from ~ 1.2 to ~ 0.7 
(IKoch et al.ll2010 ). possibly demonstrating that the role of magnetic field and turbulence evolves 
with scale. The apparent axis of the hourglass al ong the north-west south -east direction coincides 
with the axes of the bipolar outflows in CO(2-l) keto & Klaassenll2008b . CO(3-2) and HCN(4-3) 
(IShi et al.ll2010|) . Collapse and/o r accretion flow signatures, possibly orthog o nal to that, have bee n 
reported for variou s molecules in Rudolph et al.ldl990):lHo & Yound (119 96); Zha ng & Hoi (|1997|) : 
Young et all d 19981) ; Izhang et all dl998h : ISolins et all J2004l) : lKeto & Klaassenl d2008h . ~This might 



correspond roughly to the depolarization regions along the north-east south-west direction in Fig- 
ured! 

The highest resolution SMA map of W51 e2 serves as a testbed for the method developed in 
the following section. Fi gure Q] reproduces the dust continuum Stokes / map with the magnetic 
field segments (red) from iTang et al.l (|2009|) . Detected polarization segments are rotated by 90° 
in order to derive the field directions. Typical measurement uncertainties of individual position 
angles are in the range of 5° to 10°. Additionally shown in the map are features relevant for the 
method in the following section. 



3. Method 

In this section we introduce a new method which allows us to measure the (projected) mag- 
netic field strength in a molecular cloud as a function of the position in a map. Starting from an 



The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and the Academia Sinica. 
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observed correlation in Section 13.11 the details of the method and its link to observed dust con- 
tinuum Stokes / and polarization maps are given in Section 13 .21 We follow up with a qualitative 
analysis of the features of the method in Section [331 



3.1. Motivation: Magnetic Field - Intensity Gradient Correlation 

The Stokes / dust continuum conto urs and the magnetic field segments (red) from the dust 



polarization emission (ITang et al.1 120091) are shown in Figured] for W51 e2. Interestingly, field 
segments tend to be perpendicular to Stokes / emission contours, or equivalently, polarization is 
preferentially tangential to the emission contour lines. Instead of further working with intensity 
contours, we alternatively choose to investigate the intensity gradients (blue thick segments in 
Figure [TJ). The top panel in Figure |2] displays the magnetic field position angles (P.A.s) versus the 
intensity gradient P.A.s. A close correlation becomes manifest. Both P.A.s are defined counter- 
clockwise starting from north in a range of to 180°. The few data points (black empty squares) 
which originally do not seem to follow the trend, belong to pairs where both segments are close to 
P. A. = (vertical orientation), but with one of the P.A.s slightly rotated to the left and the other one 
slightly rotated to the right hand side of the vertical. Despite being closely correlated in orientation 
(small absolute differences), these pairs would mimic a bad correlation simply due to the definition 
of the P.A.s. For these cases, the original range of the magnetic field P.A. (0 to 180°) is extended 
by adding their relative difference in P.A.s to their intensity gradient P.A.S. E.g., a pair with an 
intensity gradient and magnetic field P.A. of 175° and 5° has a relative difference of 10°, and is 
therefore displayed as (175°, 185°). This re-definition does not alter the correlation in any way, but 
is merely a consequence of the ambiguity in the conventional definition of the 180° P.A. range. 

We note that the correlation is not limited to those pairs directly in the core collapsing re- 
gion, but seemingly still holds for the most northern segments where the emission contours start 
to deviate from circular shapes. The correlation - in the definition of Pearson's linear correlation 
coefficient - is 0.95. It is important to remark that, despite this tight correlation, the distribution 
of the differences (between magnetic field P.A.s and intensity gradient P.A.s) is non-Gaussian with 
an absolute mean of ~ 20°. A KS-test rejects the null hypothesis (differences following a Gaus- 
sian distribution) with a p-value of ~ 10~ 7 . As a matter of fact, the relative differences show a 
bimodal distribution (bottom panel in Figure [2]), which results from a systematic clockwise and 
counter-clockwise rotation of the magnetic field with respect to the intensity gradient. This further 
points toward a systematic origin of this correlation. All this leads to the conclusion that magnetic 
field and intensity gradient can not simply be radially aligned with the apparent misalignment just 
being the result of some measurement uncertainty. It is the focus of an on-going study where this 
correlation is further analyzed on a larger data sample in the context of the evolution of molecular 
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clouds (Koch et al. 2012; in preparation). 

In the following section, this observed correlation and the measured non-Gaussian differences 
in position angles of the magnetic field and the intensity gradient are further used to investigate the 
field strength. 



3.2. MHD and Dust Polarization 

3.2.1. Local Field Strength 

Molecular clouds are the sites of interactions of various forces. Gravity, magnetic fields, ex- 
panding shells such as from expanding HII regions or supernovae, accretion and outflows all add to 
the dynamics, and their relative importance can vary with the evolutionary stage of the star forma- 
tion site. Generally, magneto-hydrodynamics (MHD) equations - possibly for various molecules 
and dust with individual coupling terms - are needed to provide an accurate and complete descrip- 
tion of the system. In the following, aiming at making use of dust polarization measurements, a 
single MHD equation for the dust component is adopted^. Under the assumptions of negligible 
viscosity and infinite conductivity (ideal MHD case) we focus on the force equation: 

p(! + v.v)v = -v(p + g-pW +3 L(B.V)B, (!) 

where p and v are the dust density and velocity, respectively. B is the magnetic field. P is the 
hydrostatic dust pressure, cf) is the gravitational potential resulting from the total mass contained 
inside the star forming region. V denotes the gradient. As usual, the left hand side in Equation 
(0Q) describes the resulting action based on the force terms on the right hand side. These include 
the gradients of the hydrostatic pressure terms of the gas, the magnetic field and the gravitational 
potential together with the magnetic field tension term (last term on the right hand side). 

In the following, an interpretation of Equation (Q~|) is derived which can be matched to a po- 
larization observation as in Figure [Q We start by noting that the magnetic field tension term, 
~ (B • V)B, can be rewritten as: 

(B ■ V)B = B^-e SB = B^-e SB+ B*U (2) 



Simplifying and describing the molecular cloud dynamics with a single equation relies on the assumption of 
(local) collisional equilibrium between dust and other particles. Since the dust particles both couple to the magnetic 
field and also follow gravity and any other pressure gradient, other species will in a first approximation follow the 
same dynamics. 
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where a magnetic field line, B = Be SB , is directed along the unity vector e Sg with the generalized 
coordinate s B along the field force line. A change in direction, de SB /ds B , measures the magnetic 
field curvature l/R and is directed normal to e SB along the unity vector n. In an analogous way, 
(v ■ V)v can be reformulated with the velocity vector v = ve Sv directed along the unity vector e Jv . 

Generally, the gradient of the isotropic field pressure, ~ V (j^ j , in Equation CD is not nec- 
essarily directed along the field direction e SB , but can also have a component orthogonal to it. 
Therefore, when combining the Equations © and ©, a term ~ B-£^-e SB± remains which needs to 

be compared to the curvature term ~ B 2 jii. Both n and e Sfl are orthogonal to e Vfl , but they are not 
necessarily collinear. Assuming ds B± and R to be given approximately by the resolution of an ob- 
servation, the significance of this additional term can be estimated as \B^-\/ (B 2 ±) « Ml < l. 
Thus, the term can be omitted if the local change in the normal field component, AB±, is small 
compared to the total field strength. This holds for any smooth and slowly varying function of field 
strength. Equation (OQ) can then be transformed into: 

dv de Sv dP 86 1 , 1 

where the resulting directions of the gradients of pressure and gravitational potential are given with 
e Sp and e s ,, respectively. In the above formulation, generalized coordinates s v ,Sp,s^Ss along the 
directions of the unity vectors e Vv ,e V/J ,e S(/) ,e lB are used. Generally, these vectors are not collinear, 
but have different three-dimensional directions with the condition that the vector sum of the right 
hand side adds up to the resulting direction on the left hand side. Each term in Equation © is 
intentionally written with its corresponding unity vector, defining its specific direction. This will 
later be used when interpreting Equation §3§ together with an observed map. Equivalently, all 
unity vectors can be expressed in any regular coordinate system (e.g. spherical coordinates r, 9, 6) 
should this be of advantage. The unity vectors e Sv ,e sp ,e Sij) ,e SB can be interpreted as volume averaged 
directions, where the size of a volume element will depend on the resolution of an observation. In 
the above derivation we have neglected the partial time derivative, d/dt, assuming stationarity. 

Equation © describes the basic interaction of forces and the resulting motion. Linking some 
of these terms (in direction and strength) with observations will then allow us to solve for others. 
Our goal is to isolate the magnetic field strength B in the term -^B 2 ^n. Therefore, the remain- 
ing terms in Equation ([3]) need to be identified and characterized. Dust in star formation sites is 
reacting to all the force terms in Equation ©. Maps of observed dust emission are reflecting the 
overall result of gravity, pressure and magnetic field. Consequently, they represent a measure of 
the resulting motion, the left hand side in Equation ©. We make the fundamental assumption that 
a change in emission intensity is the result of the transport of matter driven by the combination of 
the above mentioned forces. Adopting this, it then follows that the gradient in emission intensity 
defines the resulting direction of motion on the left hand side in Equation ©. For the directions 
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of the gradients of the pressure and the gravitational potential we assume here, for simplicity, a 
spherically symmetrical molecular cloud, where the center can be identified from the peak emis- 
sion. This assumption is relaxed in SectionHl where the method is generalized to an arbitrary cloud 
shape. 

It is important to remark that the derivation so far is generally valid for 3 dimensions. Pro- 
jection effects in both the integrated Stokes / and the polarized emission inevitably are present in 
an observation. In the following we first proceed by identifying Equation © with an observed 
map in 2 dimensions. Projection effects are later addressed in Section l6\2l Figure [3] illustrates the 
further steps. The derivation of B is shown for a measured polarization (red) and intensity gradient 
direction (blue) in the first quadrant (with respect to the gravity center which here is supposed to 
coincide with the emission peak). For each measured pair of polarization and intensity gradient, 
its distance and direction from the gravity center defines the gravitational pull, p-^-e s ,. Generally, 
gravity, magnetic field and the intensity gradient directions are different. The schematic in Figure 
[^illustrates how the vector sum of the various terms in Equation © is constructed. Solving for the 
magnetic field term then relies on measurable angles in the orientations between polarization and 
the intensity gradient (a) and the difference between the gravity and the intensity gradient direc- 
tions (-0). Applying a sinus theorem to the closed triangle, V ^fJ^ = ^ ^ , leads to the expression 
for the magnetic field strength: 



where, in the most general case, all variables are functions of positions in a map. 

We note that - although illustrated for a simple, close to spherically symmetric case in Figure 
[3] with VP and pS/cf) aligned and pointing toward the same center - Equation © is generally valid 
for any directions of VP and pVcj). Thus, the above derivation is not restricted to a spherical 
collapsing core, but can be applied to any configuration where the various local force directions 
can be identified (see also Sections l4l and 16 .11) . 



As an important further outcome of the method, the angle factor |^ in Equation © has a 
direct physical meaning. With the magnetic field tension force term F B = ^ and the gravitational 
and pressure forces \F G +F P \ = \pV(j) + VP|, Equation © can be reformulated as: 




(4) 



3.2.2. Local Field Significance 




(5) 
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where we have introduced S B to define the field significance. This directly quantifies the local 
impact of the magnetic field force compared to the other forces involved (gravity and pressure 
gradient). Remarkably, the ratio is free of any model assumptions and simply relies on the two 
measured angles. It is purely based on the geometrical imprint by the various forces left in the field 
and intensity gradient morphologies. The angle factor, thus, also provides a model-independent 
qualitative criterion whether the magnetic field can prevent an area in a molecular cloud from 
gravitational collapse (E B > 1) or not (E B < 1). The field significance E » and its implicatio ns for 



the mass-to-flux ratio and the star formation efficiency are investigated in lKoch et al.l (|2012|) . 



ts imp i 
etalJ f 



3.3. Qualitative Analysis 

We briefly illustrate the expected features of our method when applied to a representative 
field and intensity configuration in a molecular cloud (Figure H]). A summary is given in Tabled] 
Field lines like (I) and (II) are expected in a core collapse when gravity is dragging inward the 
magnetic field which is flux-frozen to the dust particles. In the case of (I), the field line is close to 
the pole, little bent, and the gas is mostly moving along the field. Intensity gradients are mostly 
radial, pointing toward the center of gravity. Therefore, in the outer region (a), R « const, 5 small 
(sin a « 1), ip small (simp « ip « const) and thus E B « const. The field strength B will scale as 
~ r~ K , where k > describes a power-law profile resulting from density and gravitational potential 
gradient (Section 1531 ) . Small (random) changes in angles are possible due to local turbulence and 
can lead to some scatter. In the inner part (b), R might slightly decrease with the other parameters 
remaining roughly constant. This will lead to similar scalings for S B and B. A field line of type 
(II) is expected around the accretion direction. The outer field (a) is still little bent (R m const), 5 ~ 
const and small, ip small, and therefore S B m const and B ~ r~ K as in (la). In (b) the field lines start 
to bend more significantly (R decreases) and 5 increases (sin a decreases), tp is small (~ const) 
or possibly increasing if a a non-spherical central core is being formed. Thus, S B can increase 
as ~ -J— . The field strength is additionally modulated by ,/ which can lead to a steeper or 

sua ° J J y sin a r 

shallower profile than r~ K depending on the exact values of R and a. In the most inner section (c), 
the field radius is minimized, field and intensity gradient directions are close to orthogonal to each 
other (5 ~ | modulo some turbulent dispersion, a ~ 0), and ip is small. The force ratio E g can 

then possibly further increase locally. The field strength is thus mostly controlled by J and 
the power-law exponent k, which will lead to highest values here. Overall, for (I) and (II) we thus 
expect the field strength to generally increase with smaller radius toward the center. 

Additional features are seen in Figure \T\ which are possibly not directly linked to the above 
cases. Region (III) - corresponding to the NE patch around the offset (-0.3,1.8) - shows a group 
of parallel B field segments with parallel intensity gradients. In this outer region the collapse is 
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likely initiated, but does not yet reveal clear signatures, ip is changing, 5 ~ const (sin a « const), 
R « const, and S B can therefore increase or decrease. The field strength will generally be small 
here because of the low density. The relative field significance can still be important here, as 
gravity has not yet fully taken control to shape and align intensity gradient and field direction. 
Area (IV) corresponds to the most northern part in Figure [TJ This region is possibly disconnected 
from the main core, likely belonging to another smaller core being formed here. The E B values are 
dominated by large changes in ip when linked to the gravitational center of the main core. Due to 
the lack of a clear identification of a local gravity center (in which case ip would be much smaller), 
E B is probably overestimated for the two most northern segments. Possible solutions to this issue 
are addressed in Section 1431 The field strength is again small here due to the low density. 



This section focuses on the 2 remaining parameters in Equation ©, which are not yet ex- 
plicitly locally evaluated: the magnetic field radius R (curvature C = l/R) and the gravitational 
potential <p. In particular, the initial assumption of a spherically symmetrical molecular cloud from 
Section[3]is given up. It has to be stressed that the method in Section [3] derived from the Equation 
(OQ) remains unchanged, and is valid for both average or locally precise values in R and <p. 



From the derivation in Section [31 both angles a and ip, the gravitational potential and the 
pressure P are generally functions of positions (x,y) in a projected plane of sky map. The magnetic 
field line radius R (curvature C = l/R) will generally also change with positions. It is, however, 
not possible to define a curvature based on a single isolated polarization segment. For polarization 
observations with a sufficiently large number of segments (as e.g. in Figured)), a curvature resulting 
from two adjacent (or more) polarization segments can be calculated. 

Two neighboring magnetic field segments, separated by distance d, are assumed to be tangen- 
tial to a field line of radius R. The difference in their position angles, APA, is a measure of how 
much the direction of the field line has changed over d. The curvature C can then be written as 
(Figure [5]): 



With this definition of a local curvature, the calculation of the field strength B in Equation © is no 
longer limited to an average field curvature, but can be precise to any locally varying field structure. 



4. Generalizing to Arbitrary Cloud Shapes 



4.1. Local Curvature 




(6) 
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We remark that the above definition for C is not necessarily unique. In particular, determining 
which segments connect to a field line can be non-trivial. Additionally, for some neighboring 
segments, as it can also be seen in Figured! it is obvious that they can not be both tangential to 
the same field line. The curvature as defined in Equation © is therefore rather a local averaged 
curvature. In this view, a local curvature can not only be defined with two segments which are 
roughly lining up, but also with two segments side by side. The local field radius map in Figure [61 
middle left panel, is calculated by averaging the local curvatures between two adjacent segments 
over all available closest neighbors. A closest neighbor is defined through the grid spacing of the 
map. Depending on the resolution and the detection of a polarization structure, the definition of C 
can possibly be adjusted or refined. 

For W51 e2, the local field radius derived from Equation © is in the range between ~ 1" and 
~ 17" (Figure [6l middle panels). Whereas most areas in the map show rather uniform patches with 
R ~ I" - 3", regions in the north-east and south-west around the accretion plane show larger radii. 
This is possibly a consequence of the accretion flow which tends to align field lines until some 
limiting inner radius where the field will bend. Yet more likely, this is a boundary effect because 
only two almost parallel segments are available here for averaging. This leads to unphysically large 
radii. The local radius averaged over the map is < R >^ 3 "3. Discarding the extreme values gives 
<R>^ I". 



4.2. Local Gravity Direction 

Unless the gas distribution in a molecular cloud is (close to) perfectly symmetrical in az- 
imuth around an emission peak, V</> will differ in direction and strength from a simple radius-only 
dependence in the gravitational potential. Thus, for an arbitrary mass distribution, the direction 
and strength of the local gravitational pull, p§^e Scj> in Equation ©, can not any longer simply be 
constructed with the direction and distance to the reference center (emission peak in Figure Q]). 
Instead, at any position in a map, the distribution of all the surrounding mass has to be taken into 
account to calculate the resulting local gravitational force. Since the method in Section |3] is based 
on measuring angles between various force directions (including gravity) at all positions in a map, 
any substantial local deviations from a simple global gravity model can be relevant. 

For a discrete mass distribution (m,), the resulting total potential at a location r is: 0(r) = 
~GJ2'i=i irHrT' wnere G is the gravitational constant. We assume that a detected dust emission 
(like in Figured]) is proportional to the gravitating mass, and that it further traces reasonably well 
the overall gravitating potential. Additionally, it is assumed that the dust temperature is roughly 
constant over the emission map. In this case, there is no weighting for m, when linking it to the 
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integrated dust emission /,-. The local gravity force (per unit mass), g(r) = — V0(r) is then: 



where f t is a conversion factor linking dust emission and total mass, and e, is the direction between 
each detected position with dust emission (r,) and the local position (r) where the gravity force is 
evaluated. 

We remark that, if ft = / over the entire map, the local gravity direction can be derived without 
prior knowledge of /. For a typical observation, the mass (intensity) distribution is obtained by 
pixelizing a map with the beam resolution. Figure [61 bottom left panel, shows the deviations from 
an azimuthally symmetrical potential for W51 e2 when calculating the local gravity directions 
following Equation © for ft = f. This analysis shows that, although W51 e2 (Figured]) appears 
to be a rather spherically symmetrical collapsing core, the deviations are systematically positive or 
negative on one or the other side of the hourglass axis. Minimum deviations are roughly along the 
south-east north-west direction which is consistent with the hourglass axis. The average absolute 
deviation is w 9°. As further discussed in the Sections |6~4l and 1631 already this simple symmetry 
analysis points toward a systematic influence of the magnetic field. 

These deviations provide a correction for the angle ip (angle between the gravity direction 
and the intensity gradient, Figure [3]), which originally is calculated assuming an azimuthally sym- 
metrical potential (Figure [7J top left panel). Consequently, also the magnetic field strength B in 
Equation © needs to be recalculated. The angle a, as defined relative to the intensity gradient, 
remains unchanged with this local gravity correction. 

We finally note that further refinement of the local gravity correction introduced here is pos- 
sible. Local over-densities - which follow the flow of the global dynamics - can eventually locally 
collapse once they pass the Jeans mass. The north-west extension in Figure [T] is possibly such a 
case. Using Equation d7J, however, the resulting local gravity direction will always be dominated 
by the main mass (emission). Consequently, the angle ip might be over- or underestimated. This 
defect can be corrected by identifying local collapse centers with a local maxima search and a 
suitable threshold criterion. 



This section applies the method developed in Section [3] to the case of W51 e2. Starting from 
the expression for the magnetic field strength B in Equation ©, results are first presented for 
separated factors and then for the final field strength. Whereas in the original observed map in 




(7) 



5. Results 
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Figure[T]the magnetic field is gridded to half of the synthesized beam resolution, the maps in the 
following are interpolated for an enhanced visual impression. This can emphasize features along 
the field lines where polarization is detected at equally separated spacings. On the other hand, it 
might lead to non-physical values in depolarization zones along the south-west north-east direction 
or in the center, where no data are available. 



5.1. Angle Factors 

5.1.1. simp - map 

The angle ip measures the difference between the orientations of the intensity gradient and 
gravity. Small values in ip (or simp) therefore indicate that a change in the local intensity structure 
is closely following global gravity. The top panels in Figure [7] show a map of simp together with 
its trends as a function of radius and azimuth. The peak emission (Figure [B is again assumed 
as the reference center, and the gravity direction is simply taken as the radial direction from the 
center to any measured intensity location. Small values of sinip - in the range between and 0.2 
- are found for most of the magnetic field segments pointing toward the center. An X-pattern is 
revealed made up from the majority of the segments. There is a trend of larger simp values in the 
north-east and south-west where likely accretion is ongoing. The magnetic field segments in the 
far north show sin^ values close to one. Here, the intensity gradient directions are very different 
(by up to 90°) from the global gravity directions. This northern patch, therefore, seems to be 
decoupled from the main collapsing core (Section [3731 case (IV)). simp averaged over the entire 
map gives < simp >^ 0.29 (< ip >^ 20°). The average decreases to about 0.19 (< ip >^ 12°) if 
the six polarization segments in the far north are ignored. A hint of an increase in sin^ with radius 
is apparent in the top right panel in Figure |7J Again, the most northern patch dominates values 
beyond a distance of ~ 1.5". No clear structure in azimuth is revealed. 

Since ip measures a deviation from gravity, mass distributions differing from an azimuthal 
symmetry can introduce corrections to the local gravity direction (Section |4~2l ). The top left panel 
in Figure [8] displays the s'mip - map calculated from the generalized gravitational potential in 
Equation (Q. The overall patterns are still similar to the uncorrected s'mip - map, but with smaller 
values in the accretion areas. As a result, the average is reduced to < sin^ >^ 0.19 (< ip >~ 12°) 
and to about 0.09 (< ip >~ 6°) when discarding the most northern patch. 
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5.1.2. s'ma-map 

The difference in orientations between the intensity gradient and the magnetic field is denoted 
with 5 (Figure [3]). The angle a is its complement to 90°. Or, in other words, a is the difference 
in orientations between a measured (projected) polarization and the intensity gradient. Values of 
a close to its maximum of 90° (sin a ~ 1), indicate that the projected magnetic field orientation 
and the intensity gradient are closely aligned. This again reflects the initially presented correlation 
in Section 13.11 The middle panels in Figure [7] show a map of sin a with its radial and azimuthal 
dependence. Around the collapsing core, the values of sin a are between 0.85 and 1 . They decrease 
to about 0.7 in the most northern area. The intensity and magnetic field orientations are therefore 
very similar over most of the map, with an average value < sin a >^ 0.93 (< a >~ 72°). Leaving 
out the most northern part gives < sina >^ 0.95 (< a >^ 74°). Possibly there is a slight trend 
of decrease in sina with radius. No clear structures appear in azimuth, a is not affected by local 
gravity corrections as outlined in Section 1431 



5.2. Magnetic Field Force Compared to Other Forces 

' siniM _ / F B 

< sin a I local \ \ Fa+Fp \ 



The angle factor (|^) /oca/ = ( \Fg+f p \ ) = measures in a model-independent way the 
local impact of the magnetic field force (Section I3.2.2I ). The top panels in Figure [6] illustrate the 
results. Values for the relative field significance, S B , range from ~ 0.1 to ~ 1.3, with an average 
over the entire map < Eg >~ 0.33. Thus, the magnetic field can balance only about one third 
of the gravitational force. Consequently, it can not prevent the core from any further collapse. 
We stress that this result is obtained without any assumption about the mass involved. The top 
panels in Figure [6] further show that the impact of the magnetic field changes as a function of 
position. Generally, the field influence seems to become weaker toward t he center where gr avity 



likely becomes more and more dominant. This result is further explored in lKoch et al.1 (120 12|) . 



5.3. Magnetic Field Strength Maps 



Maps for the magnetic field strength B are presented in this section for two scenarios: Firstly, 
only the newly introduced force ratio S B is taken into account, leaving the other parameters in 
Equation (HJ) constant. This will illustrate the influence of S B . Secondly, the field strength is 
calculated by combining S B together with an analysis of the density profile and the resulting grav- 
itational potential. A first order estimate of the field strength, comparing the gr avitational force 
F n = and the magnetic field tension F B = is given in iTang et alj (I2009|) . Since W51 e2 



is a collapsing core, F G > F B , an upper limit of the field strength, B < 19 mG, can be derived. 



-15- 



Mr(ps 220 M Q ) refers to the gas mass enclosed within a radius 7?g(~ 1")- The radius of a mag- 
netic flux tube is 1"). The mass density p is estimated with the gas volume number density 
n# 2 (~ 2.7 x 10 7 cm" 3 ). G is the gravitational constant. The above value, B < 19 mG, will serve 
as a reference for the magnetic field strength maps presented in the following. We remark that 



Equation © is consistent with the above estimate, B = y / AttRbG^- = y/4nR B pV(l), when setting 

= 1 and neglecting the pressure gradient VP. 

In the first scenario, the B map for W51 e2 is calculated by taking into account only the 
angles tp and a. Field radius, potential and density are set constant as above. Assuming local 
changes in temperature and density to be negligible compared to gravity, the pressure gradient VP 
is omitted from Equation © in the following. From the measured ranges in ip and a (Figure [Vj 
sin?/) ps 0.1 -0.9, sin a ps 0.7- 1), it is then obvious that mostly if) will dominate the range and 
structure in a magnetic field strength map (Figure [9]). Values of the factor = £ B range from 
~ 0.1 to ~ 1.3, with an average over the entire map < S B >ps 0.33 (Figure [61 top panels). The 
resulting patterns are very similar to the ones in the sin^-map (Figure top left panel). Resulting 
variations in B are then from a few mG to ps 20 mG. The map average value is < B >ps 10 mG. 
Compared to the reference value above, the correction factor |^ generally leads to a lower field 
strength. Discarding the most northern patch - where B is likely overestimated due to the lack 
of a clear local gravity center (Section 1331) - the value becomes < B >ps 8 mG. Obviously, when 
leaving curvature, gravitational potential and density constant, the radial profile of the field strength 

1/2 

is determined by E B . As the field-to-gravity force ratio decreases toward the center, consequently 
the field strength also decreases and its largest values tend to be in the outer zones (Figure [9]). This 
result needs to be compared with the following complete analysis. 

In the second scenario, the field strength is calculated by additionally taking into account 
the profiles of density and gravitational potential. As discussed later in Section 16.21 the force 
ratio S B is little or not at all affected by projection effects. Therefore, we proceed to reconstruct 
3-dimensional radial profiles assuming spherical symmetry. For that purpose, the azimuthally 
averaged dust emission profile from Figured] is deprojected with an Abel integral. The resulting 
3-dimensional radial profile (r) is well described with a functional form f(r) ~ (r+0 1 35 )3 » whereas the 
initial projected map (R) is well fit with a profile ~ (R+ o 2 5) 2 ■ The fitting is performed on profiles with 
values binned to half of the synthesized beam resolution. Assuming the total mass of the e2 core to 
be distributed similar to its dust component, the cumulative mass m(r) ~ f' f(r')47rr ,2 dr / and the 
resulting gradient of the potential ~ ^V- can be derived. Combining all together (Equation©), the 

1 /2 

magnetic field strength profile scales as B(r) ~ (T, B (r) ■ f(r) ■ ^) . Whereas the dust emission 
map allows us to derive this functional form, its absolute normalization remains to be determined. 
Adopting < S B >ps 0.33 (Section 15721) - i.e. the observed averaged magnetic field strength can 
balance only about one third of the gravitating mass M s (ps 220 M Q ) - we normalize B so that 
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< B >=< Tj B > -(19 mG) , where the latter field strength is equivalent to 220 M & (|Tang et al 



2009). This gives < B >k, 9.4 mG. The final field strength map and its radial and azimuthal trends 
are shown in the left panels in the Figures [10] and [TT] An increase in field strength with smaller 
radius becomes apparent. A profile ~ r" 1 / 2 is overlaid for illustration. A hint of an increase in 
field strength around r ~ 1.5" possibly points toward the forming core in the northwest. This is 
found by still assuming a dominating gravity center at the emission peak. It is obvious that the 
first map (Figure |9]) can not represent the correct field strength, but merely reflects the relative 
field significance. Only by adding profiles for the density and the gravitational potential can the 
magnetic field profile correctly be derived. Density and gravitational potential only (setting Eg = 1) 
would lead to a field profile ~ r~ lJ . Taking into account the decreasing field-to-gravity force ratio 
with radius, leads to the shallower scaling B(r) ~ r" 1 / 2 . With the density profile, p ~ r~ 3 , this 
connects field and density as B(r) ~ p l l 6 . Whereas the dust emission shows little variations in 
azimuth (close to spherical symmetry), the local position-dependent measurement of E B shows 
some variations. Consequently, features in the magnetic field strength map start to be revealed, with 
possibly regions of increased strength in the accretion zones in the north-east and the south-west 
directions. In order to trace the magnetic field to the very center, higher resolution observations 
will be needed. 

We remark that the projected dust emission can also reasonably well, although with larger 
residuals, be fit with a ~ R~ l profile. The deprojected profile f(f) ~ (r+0 1 25 )2 then gives a magnetic 
field profile B{r) ~ r"° 3 and B{r) ~ p 3 / 20 (not shown). This leads to about 15% lower values in 
the field strength in the center with a similar mean field strength < B 9.8 mG. All the other 
features are left unchanged. In the above final field strength map we have neglected the variations 
in the field radius and assumed it to be constant (Figure [6] top panels). This ignores two zones 
where large field radii are likely a consequence of boundary effects (Section [47Tb . 

Finally, the influence of the local gravity correction is illustrated in the right panels in the 
Figures [10] and [TTJ Generally, features tend to be smoother. The field strengths in the accretion 
zones are less pronounced and a single center becomes apparent. The average strength is about 
7.7 mG. The radial profile of the binned values remains close to B(r) ~ r" 1 / 2 with maximum values 
around 19 mG. Structures in azimuth remain similar although with smaller variations. 

We stress that, although we are making use of the total gravitating mass M R - which sets an 
upper limit for the field strength and defines its order of magnitude - only the newly derived force 
ratio E B leads to a proper normalization and opens the possibility to derive a radial profile. This 
has become possible by solving the MHD force equation by the method introduced in Section [372] 

Errors AB in the final B field maps are estimated by propagating typical measurement un- 
certainties through Equation ©. The individual polarization P.A. uncertainty of 5° to 10° is as- 
sumed for Aa. A slightly smaller uncertainty, ~ 3° to 5°, can be expected for Aip because of 
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the averaging process in the interpolation when calculating the intensity gradient. This leads to 
AB « ±^(4irRpV4>) ■ 0.05, where the numerical factor follows from the above uncertainties com- 
bined with map averaged values for < sin a > and < simp >. For any of the above scenarios and 
< B >m 9 mG, we derive an average systematic error < AB > sys ~ ±1 mG. Errors in the profiles 
for density and gravitational potential are small after binning and azimuthal averaging. They are, 
therefore, neglected here. Besides this systematic error, there is an additional statistical uncer- 
taint y in the field orie ntation resulting from the turbulent field dispersion. Adopting 5° to 10° for 
this ( Koch et al. this contributes to the final error budget with < AB > stat ~ ±1 mG. We note 



that all uncertainties in angles are smaller than the mean value 20°) of the differences between 
magnetic field and intensity gradient PA.s (Section l3TTb . This further confirms that the correlation 
in Figure [2]can not only result from systematic and statistical uncertainties. 

In conclusion, density and gravitational potential profiles determine the order of magnitude 
of the field strength. The local features result from the angle factor These features remain 

° ° sin a 

robust with typical total errors (systematic and statistical) of AB m ±(1 -2) mG. The derived field 
strengths are also consistent with the level of magnetic fields typi cally detected by Ze eman effects 



on OH masers in compact HII regions (a few mG up to ~ 20 mG. iFish & Reid (2007)). 



6. Discussion 

6.1. Comparison, Application Regime and Limitations of the Method 

The method presented in Section [3] leads to a local magnetic field strength. It, nevertheless, 
uses strictly speaking a not exactly local but larger scale property of the field, namely its curvature 
(Figure 0). The shape of a field line is the result of various forces acting on the dust particles 
coupled to the magnetic field. These can be both long-ranging forces (like gravity) as well as 
locally acting perturbations (like turbulence). In order to solve Equation (OQ), we need to determine 
the larger scale curvature. The locally turbulent field can here act as a contaminant. This is in 
contrast to, e.g. the CF method which precisely uses the field dispersion around an average field 
direction, in combination with some turbulent velocity information, to derive a magnetic field 
strength. In principle, a curvature can be calculated from a minimum of two segments (Section 
14.11 ). but the local turbulent field can introduce some uncertainty here. Assuming an isotropic 
turbulent field component, a local curvature can be derived by averaging over several neighboring 
curvatures. As long as reasonably smooth and connected polarization patterns are detected, this 
will reveal the larger scale field curvature. This is the case for the field morphology in Figure Q] 
where a clear signature of dragged-in field lines is apparent. Higher resolution data are then needed 
to trace the field in the inner most part of the collapsing cloud where the curvature is expected to 
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change faster within a small area. 

The method further relies on identifying a noticeable effect from the various forces in Equa- 
tion (QQ) which shape the field on large scales. When assuming flux freezing, this means that 
B 2 / (87r) < P gas + p(f>+ pv 2 . In this weak field case of flux-freezing, the field lines are forced to move 
along with the gas (at least up to some scale where ambipolar diffusion sets in). As a consequence, 
gravity, intensity gradient and field orientations can be individually determined which is needed 
to calculate the angles in Equation ©. In the strong field case, B 2 /(Sn) ^> P gas . + p<p + pv 2 , where 
matter can mostly move only along field lines, the method will yield tp ~ 0, a ~ tt/2, R — > oo, 
and therefore B — > oo. (see also Section [331 ) Unless the pressure gradient VP can be identified, 
which then allows to measure ip ^0, the method can not further constrain B (see also Table [[])■ 
Geometrically, the method is based on closing a force triangle (Figure [3]). Thus, if ip ~ or a ~ 
or if both angles are zero at the same time, the triangle is no longer defined and the method fails. 

In summary, whereas the CF method is based on the field dispersion but has to isolate a mean 
large-scale field direction, the method discussed here uses the large-scale field morphology shaped 
by various forces. Therefore, the method is more directly applicable to the weak field case. A 
successful extension to the strong field case is possible if pressure and/or gravitational potential 
gradients can be localized. Generally, coherent patches with smooth changes in field and intensity 
gradient directions are of advantage as this leads to a clear identification of the field curvature 
and the acting forces. The CF method needs both dust polarization continuum and molecular 
line velocity observations. The method here uses the dust continuum Stokes / and its polarized 
emission. Whereas the CF method provides one single statistically averaged value, the approach 
here leads to a map of field strength values. 

Finally, we conclude this section by noting that the method here is not limited to dust con- 
tinuum polarization. As long as a force triangle can be identified, the method is applicable to any 
measurement which reveals the large-scale field morphology in combination with gravity and/or 
pressure forces. Notably, molecular cloud polarization maps obtained from the Gq l dreich -Kylafis 



(GK) effect, i.e. from spectral line linear polarization (|Goldreich & Kylafislll98lL 119821) . can be 
analyzed in the same way. This might open a way to study outflows which are often not ob- 
served in dust continuum due to their relatively small masses. The only complication possibly 
rises from the difficulty that it can be non-trivial to conclude whether the magnetic field is parallel 
or perpendicular to the polarization detected from the GK effect. In an even broader context, the 
method introduced here can also be applied to field morphologies detected from polarization by 
synchrotron emission or from Faraday rotation measures. The method is, therefore, not limited to 
molecular clouds, but can potentially also be adopted for, e.g. galaxy rotation measures. 
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6.2. Projection Effects 



The results in Section [5] were derived assuming the plane of sky projected geometry in Fig- 
ure [3] Here, we address the question of how an original 3-dimensional configuration affects the 
parameters in Equation © to calculate the magnetic field strength. The schematic in Figure [T3l 
illustrates the connection. We note, that like in the 2-dimensional (2D) case, gravitational pull 
(and pressure gradients), magnetic field tension and the resulting intensity gradient have to form a 
closed triangle. The projected 2D triangle can generally be deprojected with 2 different inclination 
angles: i a and Lb, corresponding to the inclinations of the deprojected orientations of the magnetic 
field tension and the gravitational pull, respectively. 

Strictly speaking, projection effects are likely to be present to some extent in all variables in 
Equation (HJ). Having deprojected the intensity emission profile and adopted an average field radius 
(Section 1531) . we limit the discussion to the remaining factor In the following, lower indices 
'3' or '2' refer to 3-dimensional or 2-dimensional projected quantities. The deprojection of the 
side &2 of the triangle (Figure [T3l) can be written as b?, = \fb\ + ^l and tanz, fc = where Tib is the 
length of the projection along the line of sight (z-direction). Analogous relations hold for the side 
a 2 - Linking a sinus theorem in the 3D triangle, = f-, to the previously used 2D expression 



where L a ,ib^ [ _ tt/2, vt/2]. It becomes obvious that for any inclination angles, the projection effects 
cancel out as long as gravitational pull (t b ) and field tension (t a ) are inclined by the same angles, 
i.e. L a = L b . Furthermore, if L a and t fe are different, projection effects are negligible as long as i a 
and L h are small. For L a and Lb within about ±40°, the projection effect is within about ±20%. The 



flat saddle around one in the center. For extreme values, i.e. inclination angles |t a |, > 75°, 
the projection effect becomes severe (f — > or / — > oo). Besides having two totally different 
inclination angles L a and Lb, we also investigate the projection effect as a function of the difference 
between the two inclination directions, A = L b -L a (Figure [141 right panel). Since the gravitational 
pull and the field tension directions are dynamically coupled, their inclination angles are likely 
to be not completely independent. The right panel in Figure [14] shows that projection effects are 
within about ±30% for a range in the differences A up to about 100° or more, up to inclinations 
of about ±45°. For larger inclinations / grows again. These situations, however, might not be of a 
concern, because at such large inclinations of the field lines, polarization emission (perpendicular 
to the field lines) might not be detectable. The analysis might, therefore, not be applicable. 



sin -02 
sin cti 



= ji, then yields: 




(8) 



left panel in Figure[[4]displays the relevant function f(L a , Lb) 



= which shows a large relatively 



In summary, for identical inclination angles projection effects cancel out. The effect is small 
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(up to ~ 30%) or even negligible for one inclination up to ~ 45° and the other one varying within 
about +50° and -50° (range in A up to about 100°). Only for extreme inclinations larger than 
about 75° projection effects can severely affect the method. Such situations might, however, not 
be detectable. It, thus, has to be emphasized that it is unlikely that severe projection effects are 
present if regular polarization patterns are observed. Consequently, the method introduced here 
leads to a total magnetic field strength. 



6.3. Shortcomings of the Method 

Without any further modeling, the method does not include components such as e.g., rotation 
or temperature variations. Effects due to rotation are not explicitly considered in Equation (Q])- 
Nevertheless, an additional term due to rotation on the right hand side of Equation (OQ) would 
primarily lead to a change in the resulting direction of motion, pv-^-e Sv + pv 2 ^- in Equation ©, 
which was identified with the intensity gradient direction. Technically, this is equivalent to an 
additional inclination in the intensity gradient and can thus be approximated as a projection effect. 
Similarly, the directions of the magnetic field lines dragged along with the rotating gas would be 
subject to an additional projection effect. As discussed in Section l6\2l the resulting uncertainty in 
the field strength is relatively small for a large parameter space of inclination angles. Alternatively, 
if the rotation velocity and axis are known, Equation (0Q) could be written in a rotating frame with 
an additional contribution from a centrifugal force. 

Density and temperature combined make up for the detected dust emission. Significant vari- 
ations in temperature could mimic mass concentrations which might then lead to wrong gravity 
directions and, therefore, to wrong magnetic field strength values. If concentrated locally such as 
e.g., in the case of radiation pressure from a star, this is unlikely to alter the global features in the 
field strength over an entire map. Similarly, small shifts in the position of a dominating emission 
peak (like in Figure [Q) introduce only minor variations in the field strength. Only if there are very 
irregular and significant changes in temperature over a map, the method is likely to fail. Additional 
information (such as from spectra) might then help to isolate the temperature effect. 



6.4. Dynamical Importance of the B Field: Flux-Freezing and the Magnetic Field - 

Intensity Gradient Relation 

Following up on Section 16.11 we expect the strong and weak field case of flux-freezing to be 
present with different characteristics. In the strong field case - gas moving along the field lines - 
the intensity gradient direction is presumably closely aligned with the magnetic field orientation. 
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On the other hand, in the weak field case - field lines dragged along and shaped by the forces 
acting on the gas - intensity gradient and field orientations can be at any angle. In the extreme 
case they can be at a 90° angle, when e.g. a gas flow is bending a field line (case (lie) in Figure 
HI). The angle 6(= n/2-a), as previously introduced in Figure |3l measures the absolute difference 
in orientations between the intensity gradient and the magnetic field. Since the cases of weak and 
strong§ flux-freezing categorize the relative importance of the field, we propose here the angle 5 as 
a diagnostic measure for the dynamical role of the B field. Small angles (~ 0) will point toward 
a dominating field, whereas larger and larger angles mean that the field is being overwhelmed by 
other forces. 

The map in Figure [[2] displays 5, revealing some areas with 5 ~ 0. These areas roughly 
coincide with the expected hourglass at larger radii. Here, the gas is possibly collapsing along the 
field lines. Additionally, it is likely that gravity has initially led to this field and intensity gradient 
configuration, with £ ^ being a mere consequence of that. Both at smaller radii toward the center 
and away from the hourglass toward the accretion plane, values for 5 increase. This might be 
an indication that the field is dynamically less important here and that gravity has largely taken 
control. It needs to be stressed that this does not necessarily mean that the field strength is smaller 
here, because B further also depends on R, V0 and VP (see also previous footnote). Moreover, 
the field can still provide pressure support which can slow down the gravitational collapse. This 
picture - with generally smaller 5 values around the core and larger values toward the north-west 
corner - is consistent with the angle factor |^ (Section [3.2.21 and top panels in Figure [6]), which 
identifies the field force to be less significant in the core. Figured!] is necessarily complementary 
to the middle panels in Figure [71 because 6 = ir/2-a. 

We have to note that the angle 5 (like a) is subject to a projection correction. In the notation 
of Section |6\21 5 3 = 7r/2-a 3 with sina 3 = |^ • sma 2 - f~ l {i a , tb)- In the particular case of a close 
to spherical collapse, |^ « 1, because both the 3-dimensional and the projected intensity show a 
clear azimuthal symmetry (FigureQ]). Assuming f(i cn L b ) to be approximately constant over a map 
(Section |6\21) , then gives sina^ « const ■ sina 2 . Thus, despite the projection effect, the dynamical 
importance of the field can still be understood from the relative differences in 5. The significance 
of the angle 5 as a diagnostic measure is further analyzed on a larger data set together with the 
correlation from Section [3TT1 in Koch et al., 2012 (in preparation). 



We remark that the terms 'weak' and 'strong' case of flux-freezing are relative to the other competing forces 
(gravity, ram pressure, etc.). They are further limited to a local region where these forces are evaluated on a map. 
Therefore, the B field in a weak field region is not necessarily weaker than the field in a strong field zone. 
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6.5. Magnetic Field and Gravity Alignment: sinw - map 

The method developed in Section [3] is based on identifying the three orientations of gravity, 
the magnetic field and the emission intensity gradient. In particular, two angles are relevant in order 
to isolate and calculate the magnetic field strength B: the angle tp between the gravity direction 
and the intensity gradient, and the angle a between the polarization orientation and the intensity 
gradient. Given the three orientations, an additional angle, u, can be read off from a map. to 
measures the difference in orientations between gravity and the magnetic field (Figure©. Thus, it 
is a measure for how effectively gravity is able to shape the field lines and aligning them with the 
direction of the gravitational pull. Small angles (small values of sinaj) indicate close alignment. 
The bottom panel in Figure [7J shows sino; for the case of W51 e2, assuming the center of gravity 
to be at the emission peak. The bottom panels in Figure [8] display again sinu; when a local gravity 
direction (Section 14.21) is adopted. The two figures are qualitatively similar. Applying a local 
gravity direction seems to sharpen the contrasts. The closest alignment (sinu; 0.1-0.2) is found 
at the outer radii close to the accretion plane in the west and in the north-west direction around the 
axis of the hourglass -like structure. An obvious poor alignment is in the northern area. Averaged 
over the entire map and excluding the northern patch we derive < s'mcu 0.21 (< u >~ 13°) 
and < s'muj >^ 0.14 (< u >~ 8°), respectively. Interestingly, after the local gravity correction, 
there is a very close resemblance between the sincj-map and the sin<5-map (Figure [12]). Since 
co measures the magnetic field - gravity alignment and 5 measures the magnetic field - intensity 
gradient alignment, the very similar maps, of course, indicate that the intensity structure is mostly 
resulting from gravity. This is also shown with the small values in the sin^-map (Figure [8l top 
panels). This is, of course, a consequence of the geometry in Figure [3] because all the 3 angles 
are dependent: to = n/2-a-ip with a = tt/2-5. Nevertheless, the clear systematic differences 
in structures between the sina; and the sin^-map (Figure [8]) lead to the conclusion that the B field 
- at least in some regions - has kept its own dynamics and is not yet aligned and overwhelmed 
by gravity. This result is then also consistent with the finding of zones of weak and strong flux- 
freezing (Section l6\4l) which characterizes the dynamical role of the magnetic field. The systematic 
influence of the field is further manifest in the systematic deviations from an azimuthally symmetric 
emission structure, seen in the bottom panels in Figure[6] Projection effects (Section I6T21 are again 
likely to change only the absolute numbers, leaving relevant features and comparisons between 
maps still intact. 

7. Summary and Conclusion 

Dust polarization observations in molecular clouds often show position angles of magnetic 
field segments which are close to perpendicular to the dust continuum intensity contours. This 



correlation does not seem to be arbitrary, but can be interpreted in the context of ideal magneto- 
hydrodynamics (MHD). Based on this, we put forward a new method to derive a local magnetic 
field strength. The key points are summarized in the following. 



1. New method for magnetic field strength mapping: A new method is proposed which leads 
to a position-dependent estimate of the magnetic field strength. The approach is based on 
measuring the angle between the polarization and the Stokes / intensity gradient orientations. 
Assuming that the intensity gradient is a measure for the resulting direction of motion in the 
MHD force equation, this angle quantifies how tightly the direction of motion is coupled to 
the field orientation. In combination with the angle in between the orientations of gravity 
(and/or the pressure gradient) and the intensity gradient, the two angles can be linked to the 
field strength by identifying a force triangle. 

2. Application regime: The method relies on separating noticeable effects from various com- 
ponents in the MHD force equation, in particular from gravitational pull and field tension. 
Consequently, the method is more directly applicable to the case of weak field flux freezing, 
where the field morphology is clearly shaped by these forces. Collapsing (spherical) cores 
in star formation sites are of most immediate interest because the directions of gravity and 
magnetic field are relatively easily determined. Nevertheless, the technique is generally valid 
and applicable to any measurement (also strong field flux freezing), where the directions of 
magnetic field, pressure gradients and/or gravitational pull can be localized. With the con- 
cepts of local curvature and local gravity direction, the method can be further generalized. 
Moreover, the method is not limited to dust polarization measurements, but can potentially 
also be applied to, e.g. molecular line linear polarization maps from the GK effect or even 
to Faraday rotation measure maps for galaxies. 

Here, the method is applied to SMA dust polarization data of the collapsing core W51 e2. 
Averaging over the entire core, a field strength of ~ 7.7 mG is found. Variations in the field 
strength are detected with an azimuthally averaged radial profile B(r) ~ r" 1 / 2 . Maximum 
values toward the center are about 19 mG. The current data lack resolution in order to probe 
the innermost part of the core, where the largest field strength is expected. 

3. Comparison with Chandrasekhar-Fermi ( CF) method: The method introduced here uses the 
large-scale field curvature. The local turbulent field dispersion is acting as a contaminant. 
This is in contrast to the CF method which precisely uses the field dispersion - after separat- 
ing a mean field orientation - in combination with a turbulent velocity information in order 
to derive a field strength. Whereas the CF method yields one statistically averaged field 
strength for an ensemble of polarization segments, the method here provides a field strength 
at each location of a detected polarization segment. 
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4. Limitations and shortcomings of the method: Technically, the method fails where the force 
triangle can not be closed. This leads to the unphysical value B — > oo. A triangle degenerat- 
ing to a single line corresponds to the situation where the field is very weak and, therefore, 
the intensity gradient does not deviate from the gravity direction. The method then still 
consistently yields the limit B — > 0. 

The method is not entirely free of projection effects. In the most general case, the force 
triangle has to be deprojected with two different inclination angles. Nevertheless, the projec- 
tion correction enters as the ratio of the two inclination angles. Thus, in the most favorable 
case with identical inclination angles for gravity and field orientations, the projection effect 
cancels out and the method leads to a total magnetic field strength. Even if the inclination 
angles differ from each other, the possible error due to unknown inclinations is relatively 
modest, within ±30% to ±50%, for a large parameter space. In any case, if the inclina- 
tion angles do not vary significantly within the same source, the method always provides a 
way to detect relative differences in the field strength. The accuracy and reliability of the 
method are limited by how well projection effects and forces defining the field morphology 
can be identified and quantified. Additional effects like, e.g. rotation and local variations in 
temperature are likely of minor importance. 

5. Weak and strong flux freezing: The dynamical role of the magnetic field - in the sense of 
weak or strong flux freezing - can be assessed through the angle between the magnetic field 
and the intensity gradient orientations. Angles close to zero point toward a dominating role 
of the field (or an overall controlling gravitational pull). Larger angles indicate that the field 
lines are more significantly shaped by other forces. Although this criterion alone does not 
constrain the field strength, it provides a measure to characterize the role of the magnetic 
field as a function of position in a map. 

6. Dynamical significance of magnetic field: The method additionally provides a model-independent 
criterion to quantify the magnetic field force compared to the sum of any other forces in- 
volved. As a result of the force triangle, the local significance of the magnetic field can be 
evaluated simply with the ratio of the angle between intensity gradient and gravity direction, 
and the angle between intensity gradient and polarization orientation. These angles reflect 

the imprint of the various forces onto the dust and polarization morphologies. This measure 
can then constrain the local impact of the magnetic field and quantify its support against 
gravitational collapse. 

The authors thank the referee for valuable comments which led to further significant insight 
in this work. P.T.P.H. is supported by NSC grant NSC97-21 12-M-001-007-MY3. 
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Fig. 1. — The collapsing core W51 e2 (jTang et al.l 120091 ) . Contours show the Stokes / dust con- 
tinuum at 0.87 mm with a synthesized beam resolution of about O'.'l (gray ellipse in the lower 
left corner) with natural weighting. Colors correspond to the color wedge on the right hand side 
with units in Jy/beam. Contours are linearly spaced in steps of 130 mJy/beam. Overlaid are the 
magnetic field segments (thick red segments) at the locations where polarized emission was de- 
tected. Magnetic field segments are plotted by rotating the polarization segments by 90°. Only 
polarization data with a flux above 3a lp , the rms noise of the polarized intensity, are included. 
The blue vector field displays the gradient directions of the dust continuum emission, with most 
vectors pointing toward the emission peak. Highlighted are those emission intensity gradient di- 
rections (thick blue segments) at the locations of the magnetic field segments. The length of the 
segments is arbitrary, and for visual guidance only. The axes offset positions (in arcsec) are with 
respect to the original phase center of the observation at Right Ascension (J2000)= 19 h 23 m 43 s .95, 
Declination (J2000)=14°30'34."00 
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Fig. 2. — Top panel: the correlation (blue filled squares) between the intensity gradient PA.s (blue segments in 
FigureHJ and the magnetic field P.A.s (red segments in Figure [T). The black empty squares belong to pairs with P.A.s 
close to P.A. = 0, one P.A. being slightly rotated to the left and the other one to the right hand side of the vertical. In 
order to properly display their correlations, the magnetic field P.A. is re-defined beyond the to 180° range for these 
cases (blue filled squares above 180° and below 0°). For visual guidance added is the straight blue line, representing 
a perfect correlation. Also shown are the absolute differences (< 90°) between the PA.s for each pair (red crosses). 
Both PA.s are defined counter-clockwise starting from north. Bottom panel: the bimodal distribution of the relative 
differences (in the range between -90° to +90°) between the intensity gradient and the magnetic field PA.s, reflecting 
a systematic clockwise and counter-clockwise rotation of the field with respect to the intensity gradient. 
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Fig . 3 . — Illustration of the Equation with the relevant angles in order to solve for the magnetic field strength B. A 
pair of magnetic field (red solid line) and emission intensity gradient segments (blue solid line) are shown at a location 
( Ax, Ay) with respect to a common reference center. The originally measured polarization orientation (rotated by 90° 
with respect to the field orientation) is indicated as a red dashed line. Position angles are measured counter-clockwise 
with respect to north (N). S < it/2 measures the angle in between the intensity gradient and the field orientation, a is 
its complement to ir/2. (W+ pV<f>) are assumed to be known in direction and strength. For simplicity, as suggested 
from Figure [T] for the case of W51 e2, the emission peak is assumed as the gravity center, and thus defines direction 
and strength. Here, it also serves as the reference center, ip denotes the deviation of the intensity gradient from the 
direction of {VP+ pV<j>). The corresponding deviation for the magnetic field is measured with the angle to. The force 
triangle - (VP+ pV0) pointing toward the gravitating center, the resulting intensity gradient orientation (left hand side 
in Equation (f3]l) and the direction of the restoring magnetic field tension (normal to the field segment) - can then be 
closed by intercepting the red dotted and black dotted lines, where the former one is shifted parallel to the red dashed 
line. 
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Fig. 4. — Qualitative analysis (Section [3731) : Illustration of the expected features of the method. 
The black solid lines are intensity emission contours. Blue and red segments are intensity gradient 
and magnetic field orientations, respectively, as in the Figures [T] and [3] The labeling refers to the 
cases discussed in Section [331 and listed in Table [TJ All cases have similar equivalents in the W51 
e2 map in Figured] 



Fig. 5. — Illustration of the local curvature derivation. Position angles (PA.s), separated by an ob- 
served distance d, are measured counter-clockwise with respect to north (N). The 2 magnetic field 
segments (solid red lines) are assumed to be tangential to a connecting field line (thick solid black 
line) with radius R. The angle £ measures the difference in between 2 PA.s, £ = \PA.\ -PA. 2 \ = 
APA.. R is read from the isosceles triangle with COS7 = where 7 = = |(7r- APA.). This 
then leads to Equation ©. Note that the derivation uses the measured relative quantities d and 
APA.. Technically, this allows to detach the magnetic field segments from their absolute posi- 
tions, by keeping their relative separation and APA., and ensures that a connecting field line can 
be constructed. Strictly speaking, even this definition of curvature is not exactly local because it 
involves at least 2 segments separated by a distance d. 
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Fig. 6. Position dependence of parameters in Equation {4). Magnetic field segments are shown in white. Top left panel: the combined factor 

= Eg, which quantifies the significance of the magnetic field force compared to the other forces (Section |3.2.2) . ip is measured with respect 
to a spherically symmetrical gravitational potential. Middle left panel: the local magnetic field radius calculated from the method in Section |4~T1 
Units are in arcsec. Bottom left panel: the difference in the direction of the gravitational pull between a spherically symmetrical potential and a 
local gravity direction as defined in Section |4.2| Units are in degree. The difference is shown relative to the symmetrical assumption, revealing a 
symmetry axis along the southeast- northwest direction. The right panels illustrate the tendencies for Eg, B field radius and difference in gravity 
direction (from top to bottom) as a function of radius (distance from emission peak) and azimuth angle (measured counter-clockwise from west). 
Mean values averaged over the entire maps are shown with the red solid lines. For the local gravity difference the mean is calculated from the 
absolute differences. 
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Fig. 7. — Relevant angles for the derivation of the magnetic field strength, calculated for W51 e2 (Figure [TJ. 
Magnetic field segments are shown in white. Top left panel: sin^-map, where rp (< ir/2) (Figure [3) is the difference 
between the orientations of the intensity gradient and gravity. Middle left panel: sina-map. Tr/2-a (< n/2) measures 
the deviation between the magnetic field and the intensity gradient orientations. Bottom left panel: sinw-map, where 
uj is the difference between the magnetic field and gravity orientations. The right panels show the corresponding 
tendencies for each angle as a function of radius (distance from emission peak) and azimuth angle (measured counter- 
clockwise from west). Mean values averaged over the entire maps are shown with the red solid lines. 
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Fig. 8. — Angles corrected for local gravity directions. Magnetic field segments are shown in white. Top panels: 
sin^ measuring the deviation between the local gravity and the intensity gradient orientations. Bottom panels: sinw 
measuring the deviation between the local gravity and the magnetic field orientations. 
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Fig. 9. — The magnetic field strength (B) derived by taking into account only variations in the 
field-to-gravity force factor |^ = S B . Units are in mG. All the other parameters in Equation ([4]) 
are constant over the map. Magnetic field segments are shown in white. The right panels show the 
field strength tendencies as a function of radius (distance from emission peak) and azimuth angle 
(measured counter-clockwise from west). Mean values averaged over the entire map are shown 
with the red solid lines. 
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Fig. 10. — Final magnetic field strength maps calculated from Equation (HJ). Units are in mG. 
Taken into account are the deprojected profile for the density, the resulting gravitational potential 
gradient profile and variations in the force ratio S B (left panel). In the right panel, corrections due 
to the local gravity direction are additionally considered. Mean field strengths are about 9.4 mG 
and 7.7 mG for the left and right map, respectively. The field radius is set constant. Magnetic field 
segments are shown in white. Typical measurement uncertainties lead to an average systematic 
error in the field strength of < AB > sys ~ ±1 mG. 
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Fig. 11. — Radial and azimuthal trends for the magnetic field strength maps presented in Figure 
ITOl Radial profiles ~ r" 1 / 2 are overlaid for illustration in the top panels. Mean values averaged 
over the entire map are shown with the red solid lines. 
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Fig. 12. — The angle 5 (absolute difference in between the intensity gradient and the magnetic 
field orientations (Figure [3])) as a diagnostic tool to distinguish between weak and strong case of 
flux-freezing. Magnetic field segments are shown in white. The values for sin 6 are displayed in 
the left panel in order to facilitate a comparison with the sin a (Figure 13 middle panels) and the 
s'mu) maps (Figured bottom panel; Figured bottom panel). 
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Fig. 13. — Projected and deprojected triangles for the method presented in Section[3] Lower indices '3' or '2' 
refer to 3-dimensional or 2-dimensional projected quantities. The projected triangle (dashed lines) in the xy-plane is 
deprojected along Hb and ir a (parallel to the line of sight) to form a 3D triangle (solid lines). The angles ip2 and ai are 
identical to t/j and a in Figure [5] Lb and i a are the inclination angles of the gravitational pull and the magnetic field 
tension, respectively. 
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Fig. 14. — Left Panel: The projection function f(i a ^b) = ^377^ for the two inclination angles i a 
and Lb (Figure [13]). The relatively flat saddle in the center shows that projection effects are small 
(f ~ 1) for L a and tb within about ±40°. Highlighted in black and red are the parameter spaces 
where the projection effects are within ±30% (f within 1 ± 0.3) and within ±50% (f within 
1 ±0.5), respectively. For much larger inclination angles polarization is not easily detectable. It is, 
thus, unlikely that severe projection effects are present if regular polarization patterns are observed. 
Right Panel: The projection effect as a function of the difference A between the two inclination 
directions, A = L h -i a , for various fixed inclinations i a . The dashed line marks / = 1. A ±30% 
bound is shown with the dotted lines. 



Table 1. Qualitative Analysis 



Eg B Comments 



General 



4-t 

¥0 



¥0 

if 



it it, ~ ^simp-R ■ r 
ti ti, ~ \f^-r- K 



Eg and B 4" with larger angle rjj 
Eg and B f with smaller angle o 



Specific Cases 



(I) (a) 


small, ft: const 


7T 

~ 2 


ft; const 


ft; const 


ft; const- r -K 


field line close to pole, Eg ft: const, B f with smaller r 


(b) 


small, RJ const 


7r 

~ 2 


ft; const 


ft; const 


ft; const- r _K 




(II) (a) 


small, ft; const 


7T 

~ 2 


ft; const 


ft; const 


ft; const- r~ re 


field line around accretion direction, Eg ft; const, B f with smaller r 


(b) 


small, ft const 




4- 


~ 1 

sin a 


\^ sin a ^ 


Eg f , 6 f with smaller r steeper or shallower than (depending on R and a) 


(c) 


small, ft; const 


4-, ~o 


4-' Rmin 


~ 1 

sin a 


sin a ^ 


Eg "f, maximum field strength B 


(HI) 


n 


ft: const 


ft: const 


ti 


t4-> ~ v /sm V' ' r " 


outer area, no clear signatures yet, B generally small due to r~ K 


(IV) 


large 


ft; const 


ft: const 


large 


~ V sini/) • r~ K 


disconnected area, no clear gravity center localized, B small due to r~ K 


Limits 
















^0 


/0 


(i)R 1 , (ii) -> oo 2 


^0 


(i) -> 1 , (ii) -> oo 2 


triangle degenerates to single line 




¥0 


-> 


(i) R, (ii) -> oo 


— > oo 


(i) — > oo, (ii) — > oo 


method fails, triangle not closed 




-> 


-> 


(i) R, (ii) -> oo 


-M 


(i) ~ v^, (ii) -> oo 


method fails, triangle not closed and degenerated to single line 



Note. — Qualitative analysis of the force ratio Eg = and the magnetic field strength B= y jjS|^ pV</> ■ 4ttR (Section [3.3t . The pressure gradient is neglected here. The profile of the 
density and the resulting gradient of the gravitational potential is represented with r~ K with k > (Section [5.3l . The specific cases labeled (I) to (IV) refer to Figure|4] The cases (lb), (lib) 
and (lie) are not yet observed due to the lack of higher resolution data. Down-arrows (4-) and up-arrows (y) indicate decreasing and increasing values, respectively, 'ftconst' indicates that 
small local variations can occur which lead to some scatter in Eg and B (Figure [TO . Failure of the method - triangle can not be closed - refers to Figure|5]. 

'The case if) — > states generally that the intensity gradient does not deviate from the gravity direction for any orientation of the field lines. Therefore, the field strength must be insignificant, 
i.e. B — > 0, which is also consistently reproduced by the method. 

2 The case a ¥0 ~ n/2 with R — > oo - i.e., straight B field lines aligned with the intensity gradient - represents the strong field flux freezing situation (see Section l6Tt and consistently 
yields a dominating field strength (B — > oo, Eg — > oo). 



